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Abstract 

We consider the Cartesian product X of n finite intervals of integers and a map F 
from X to itself. As main result, we establish an upper bound on the number of fixed 
f/j ' points for F which only depends on X and on the topology of the positive circuits 

^, ' of the interaction graph associated with F. The proof uses and strongly generalizes 

a theorem of Richard and Comet which corresponds to a discrete version of the 
^^ \ Thomas' conjecture: if the interaction graph associated with F has no positive 

Q^ ' circuit, then F has at most one fixed point. The obtained upper bound on the 

CN ■ number of fixed points also strongly generalizes the one established by Aracena et 

^~~r ' al for a particular class of Boolean networks. 

^^ \ Key words: Discrete dynamical system. Discrete Jacobian matrix. Interaction 

^— V \ graph. Positive circuit. Fixed point. 
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. 5t , 1 Introduction 



We are interested by tlie number of fixed points for maps tliat operate on the 
Cartesian product of n finite intervals of integers (when this product is {0, l}", 
such maps are usually called Boolean network). Our motivation comes from 
biology, where these maps are extensively used to describe the behavior of 
gene networks. The context is then the following. 

When studying gene networks, biologists often illustrate their results by inter- 
action graphs. These are directed graphs where vertices correspond to genes 
and where edges are labelled with a sign: a positive (resp. negative) edge from 
i to i means that the protein encoded by gene j activates (resp. represses) the 
synthesis of the protein encoded by gene i. These graphs are then used as basis 
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to generate dynamical models describing the temporal evolution of the concen- 
tration of the encoded proteins (see [1] for a literature review). Unfortunately, 
these models require, in most cases, unavailable information on the strength of 
the interactions. One is thus faced with the following difficult question: Which 
dynamical properties of a gene network can be inferred from its interaction 
graph (in the absence of information on the strength of the interactions)? 

In this paper, we focus on this question in a general discrete modeling frame- 
work. The set of states of a network of n genes is represented by the product 
X = nr=i -^i of n finite intervals of integers. Each interval Xj then corresponds 
to the set of possible concentration levels for the protein encoded by gene i. 
On one hand, the dynamics of the network is described by the successive iter- 
ations of a map F from X to itself whose fixed points correspond to the stable 
states of the network. At this stage, it is worth noting that the number of 
stable states is a key feature of gene networks dynamics: according to an idea 
of Delbriick |2j, the presence of multiple stable states is one possible mecha- 
nism for biological differentiation. One the other hand, the interaction graph 
of the network is deduced from F in two steps. First, to each state x E X 
and to each directional vector v G {—1,1}" such that x + f G X is associated 
a local interaction graph Gp{x,v) which contains a positive (resp. negative) 
edge from j to i if 

Ji[Xi, . . . , Xj -\- Vj, . . . , Xji) Ji\^) 

is positive (resp. negative) (/j denotes the ith component of F). Then, the 
global interaction graph G{F) of the network is defined to be the union of 
all the local interaction graphs. Note that each local interaction graph is a 
subgraph of the global one, and that the global interaction graph can have 
both a positive and a negative edge from one vertex to another. 

In this setting, Richard and Comet P| partialy answer the previous question 
by proving a well known conjecture of Rene Thomas relating the stable states 
of the network to the positive circuits of its local interaction graphs (a circuit 
is positive if it has an even number of negative edges). A weak form of their 
result follows (the original statement needs additional definitions and is given 
latter in the paper): 

Theorem 1 [|3j Let X be a product of n finite interval of integers, and let 
F be a map from X to itself. If all the local interaction graphs Gf{x,v) are 
without positive circuit, then F has at most one fixed point. 

Aracena, Demongeot and Goles [IfS] proved another theorem relating stable 
states to positive circuits. They establish, in a boolean context, an upper 
bound on the number of stable states which only depends on the positive 
circuits of the global interaction graph of the network: 



Theorem 2 [4ll5] Let F be a map from {0, 1}" to itself such that G{F) has 
no both a positive and a negative edge from one vertex to another. If I is a 
subset of {1, . . . ,n} such that each positive circuit of G{F) has at least one 
vertex in I, then the number of fixed points for F is less than or equal to 2'^'. 

The main resuh of this paper is a significative generahzation of both Theoreni[T] 
and Theorem [2l A weak form of this resuh is: 

Theorem 3 Let X = nr=i -^i ^^ '^ product of n finite interval of integers, and 
let F be a map from X to itself. If I is a subset of {1, . . . , n} such that each 
positive circuit of each local interaction graph Gf{x,v) has at least one vertex 
in I, then the number of fixed points for F is less than or equal to Yiiei \^i\- 

Theorem [3] imphes Theorem^ since if aU the local interaction graphs Gf{x, v) 
are without positive circuit then 7 = satisfies the conditions of Theorem [31 
So the corresponding bound is 1 and Theorem [1] is recovered. Theorem [3] also 
implies Theorem [2l Indeed, let F be a map from X to itself and suppose / to 
be such that each positive circuit of G{F) has at least one vertex in /. Then 
I satisfies the conditions of Theorem [3] because each local interaction graphs 
Gf{x, v) is a subgraph of G{F). So the corresponding bound is Ilig/ \Xi\ and it 
equals 2^^^ in the particular case where X is the n-cube {0, 1}". We thus recover 
the conclusion of Theorem [2] (even if G{F) has both a positive and a negative 
edge from one vertex to another). The proof of Theorem [31 which is done by 
induction on / with Theorem [H as base case, is independent of the proof of 
Theorem [21 given in [HI5] . Note also that Theorem [21 does not imply Theorem [H 
even if this latter is stated for maps F from {0, 1}" to itself such that G{F) 
has no both a positive and a negative edge from one vertex to another. 

The paper is organized as follow. In Section [21 in order to obtain a bound 
stronger than the one mentioned above and more relevant from a biological 
point of view, we focus on the asynchronous iterations of F that Thomas use 
to describe the dynamics of gene networks piTfSllQ] . First, we represent these 
iterations under the form of a directed graph T{F) on X usually called asyn- 
chronous state transition graph. Then, we define the attractors of T{F) to be 
the smallest subsets of states without output edges in T{F). The fixed points 
of F then correspond to particular attractors. In Section [31 we characterize a 
subgraph Gf{x, v) of Gf{x, v) which only depends on T{F) and which is, for 
this reason, well suited to the study of T{F). In Section [H we establish our 
main result: an upper bound on the number of attractors in T{F) which only 
depends on the map Of and which has Theorem [31 as immediate consequence. 
Final comments are given in Section [51 These are about the influence of con- 
nections between positive circuits and the interest of the established bound 
in the context of the so called Thomas ' logical method piTfSfQ] which is, in 
practice, one of the most usual discrete modeling method of gene networks. 



2 Asynchronous state transition graph and attractors 



Let X = IliLi Xi be the product of n finite intervals of integers of cardinality 
strictly greater than 1, and consider a map F from X to itself, 

X = (xi, ...,Xn)eX h-> F{x) = (/i(a;), . . . , /„(x)) G X. 



In the following definition, we attach to F a directed graph on X called 
asynchronous state transition graph. According to Thomas p[7f8|9] . this state 
graph can be seen as a model for the dynamics of a network of n genes: the set 
of vertices X is the set of possible states for the network (each interval Xi cor- 
responds to the possible concentration of the protein encoded by gene i) , and 
each path corresponds to a possible evolution of the system. [Asynchronous 
state transition graphs can also be seen as discretizations of piecewise-linear 
differential systems, see [XOfTT] for instance.] 

Definition 1 The asynchronous state transition graph of F is the directed 
graph T{F) whose set of vertices is X and which contains an edge from x to y 
if there exists i G {1, . . . , n} such that 

fi{x)j^Xi and y = x + sigia{fi{x) - Xi) ■ Ci, 

where Ci denotes the n-tuple whose ith component is 1 and whose other com- 
ponents are 0, and where sign(a) = a/\a\ for all integer a 7^ 0. 

[Following this description of the dynamics, fi{x) can be seen as the value 
toward which the concentration Xj of the protein encoded by gene i evolves: 
at state z, there exists a state transition allowing the ith component of the 
system to increase (resp. decrease) if and only if Xi < fi{x) (resp. Xi > fi{x)).] 

The fixed points of F have no successor in T{F) and naturally correspond to 
the stable states of the system. In the next definition, we introduce the notion 
of attractor which extends, in a natural way, the notion of stable state. 

Definition 2 A trap domain ofr{F) is a non-empty subset A of X such that, 
for all edges {x,y) ofr{F), if x & A then y & A. An attractor ofr{F) is a 
smallest trap domain with respect to the inclusion relation. 

In other words, the attractors of T{F) are the smallest set of states that we 
cannot leave. They extend the notion of stable state in the sense that x is a 
fixed points of F if and only if {x} is an attractor of T{F). Note also that 
there always exists at least one attractor (since A is a trap domain). Other 
easy observations follow: (1) From each state, there is a path which leads to an 
attractor (this is why one can say that attractors perform, in weak sense, an 



attraction); (2) Attractors are strongly connected components; (3) Attractors 
are mutually disjointed (this point used in the proof of our main result). 



3 Discrete Jacobian matrix and interaction graph 



In this section, we introduce a notion of local interaction graph well suited to 
the study of r(F). We proceed as in [3] by first introducing a discrete Jacobian 
matrix for F based on a notion of discrete directional derivative. 

Let X' be the set of couples (x,f ) such that x G X, i; G { — 1, 1}" and x+v G X. 

Definition 3 For all (x, v) G X' , we call Jacobian matrix of F evaluated at x 
along the directional vector v the nx n matrix F'{x, v) = {fij{x, v)) defined by 

fij{x,v) = '-^ (2,j = l,...,n). 



[If Vj is positive (resp. negative), then fij{x, v) may be seen as the right (resp. 
left) partial derivative of fi with respect to the jth variable evaluated at x. In 
both cases, fij{x,v) is a natural discrete analogue of {dfi/dxj){x).] 

An interaction graph is here a directed graph whose set of vertices is {1, . . . , n} 
and where each edge is provided with a sign. More formally, each edge is 
characterized by a triple {j,s,i) where j (resp. i) is the initial (resp. final) 
vertex and where s G { — 1,1} is the sign of the edge. The set of edges of an 
interaction graph Q is denoted E{Q). An interaction graph ^ is a subgraph of 
an interaction graph Q' if E{Q) C E{Q'). 

Definition 4 We call interaction graph of F evaluated at (x, v) G X' , and 
we denoted by Gf{x,v), the interaction graph which contains a positive (resp. 
negative) edge from j to i if fij{x,v) is positive (resp. negative). 

[To illustrate this definition, assume that fij{x,v) is positive and that Vj = 1. 
Then, fi{x) < fi{x + Cj) so we can say that, at state x, an increase of Xj 
induces an increase of fi, that is, an increase of the value toward which the ith 
component of the system evolves. In other words, j acts as an activator of i, 
and we have a positive edge from j to i in Gf{x, v).] 

In our context, the obvious fact that Gf{x,v) does not only depend on r{F) 
is not satisfactory since it is commonly accepted that the interaction graph of 
a network only depends on its dynamics, which is here characterized by T{F). 
This lead us, as in [3], to slightly modify the definition of Gf{x, v) in order to 
obtain an interaction graph Gf{x,v) which only depends r(F). 



Definition 5 We call interaction grapli of F evaluated at (x, v) G X' witli 
tliresliolds, and we denote by Gf{x,v), the interaction graph which contains a 
positive (resp. negative) edge from j to i if fij{x,v) is positive (resp. negative) 
and if fi{x) and fi{x + VjCj) are on both sides of (the threshold) Xi + Vi/2. 

[a and h are on both sides of c if a < c < 6 or 6 < c < a.] 

Remark 1 Gf{x,v) is a subgrapli of Gf{x,v) (often strict since the addi- 
tional condition "on both sides of the threshold" is rather strong). 

Remark 2 The introduction of Gf{x,v) has been motivated by arguments 
coming from the modeling context. Another relevant argument is the following: 
because Gf{x,v) is a subgrapli of Gf{x,v), all the incoming results remains 
valid but becomes less strong when stated with Gf{x,v) instead of Gf{x,v). 

Remark 3 In the Boolean case, i.e. when X = {0, 1}", Gf{x,v) = Gf{x,v). 

Definition 6 We call global interaction graph of F, and we denote by G{F), 
the interaction graph whose set of edges is [j{x.v)(^x' E{Gf{x,v)). 

Obviously, G{F) only depends on T{F) and can thus be seen as the global 
interaction graph of the network of dynamics T{F). Note that G{F) can have 
both a positive and a negative edge from one vertex to another. 

Now, we recall the notion of positive circuit and the notion of positive feedback 
vertex set. This has been introduced by Aracena et al [IfS] to study the fixed 
points of Boolean networks. 

Definition 7 A positive circuit in an interaction graph Q is a non-empty 
sequence of edges, say 

(jl, Si,ii), (j2, S2, «2), • • • , {jr, Sr, v), 

such that: ik = jk+i for 1 < k < r (the sequence is a path); ir = ji (the path 
is a circuit); the vertices jk are mutually distinct (the circuit is elementary); 
the product of the signs Sk is positive (even number of negative edges). 

Definition 8 [3] yl positive feedback vertex set of an interaction graph Q is 
a subset / C {1, . . . ,n} such that each positive circuit of Q has a vertex in J. 

One can remark that: (1) The set of vertices of Q is always a positive feedback 
vertex set of Q] (2) The empty set is a positive feedback vertex set of Q if 
and only if Q has no positive circuit; (3) If Q' is a subgraph of Q then all the 
positive feedback vertex sets of Q are positive feedback vertex sets of Q' . 



4 Positive circuits and attractors 



As previously, let X = OiLi -^i be the product of n finite intervals of integers 
of cardinality strictly greater than 1, and let F be a map from X to itself. 

We are interested by the relations between the map Gp (defined on X') and 
the number of attractors in T{F). The following theorem, presented in [3J as 
solution of a discrete version of the Thomas' conjecture, gives such a relation. 

Theorem 4 [3] If Gf{x,v) has no positive circuit for all {x,v) G X' , then 
T{F) has a unique attractor. 

The following theorem extends the previous one by providing, without any 
condition on the map Gp, an upper bound on the number of attractors in 
r(F) which only depends on Gp. 

Theorem 5 (main result) For each i G {1, . . . ,n}, let T^lGp) be the set of 
real numbers t for which there exists {x,v) G X' such that t = Xi + Vi/2 and 
such that i belongs to a positive circuit of Gf{x,v). Suppose I to be, for all 
{x,v) G X' , a positive feedback vertex set of Gf{x,v). Then, the number of 
attractors in T{F) is less than 



n [\T^iGI 



iei 



Proof — We reason by induction on /. Suppose / to be, for any (x, v) G X', a 
positive feedback vertex set of Gf{x,v). 

Base case. If J = it means that there is no {x, v) G X' such that Gf{x, v) has 
a positive circuit. So, following Theorem HJ T{F) has at most one attractors 
and the theorem holds. 

Induction step. Suppose that / 7^ 0. The induction hypothesis is the following: 

Induction hypothesis: Let F be a map from X to itself. If / is, for all {x, v) G 
X', a positive feedback vertex set of Gp{x,v), and if / is strictly included 
in J, then T{F) has at most Yii^f \Ti{Gp)\ + 1 attractors. 

Without loss of generality, suppose that 1 G /. Let P be the partition of 
Xi whose elements Y are the maximal intervals of Xi (with respect to the 
inclusion relation) verifying 

Vt G Ti{GF), t < min(r) or max(r) < t. (1) 



Remark that, by definition, 

|P| = |Ti(G'^)| + l. (2) 

Let Y be any interval of P, and consider the map F = (/i, . . . , /„) : X —>■ X 
defined by fi = fi ioi i > 1 and by 



Wx e X, A(x) 



min(y) if fi{x) < min(y) 
/i(x) if/i(x)GF 
max(y) if /i(x) > max(y). 



Then, for all x,y E X, 

Ux)<h{y) ^ h{x) < Ux) < Uy) < h{y) (z = l,...,n). (3) 

Indeed, this is obvious for i > 1, and for i = 1 it is sufficient to remark that 

h{x) < h{y) => h{x) < max(r) ^ h{x) < fi{x), 

and that 

Mx) < Uv) => min(r) < f\{y) => f\{y) < h{y)- 

Now, we prove that, for all {x^v) G X', 

Gp{x, v) is a subgraph of Gf{x, v). (4) 

Let (x, v) G X' and suppose (j, s, i) to be an edge oiGp{x, v). According to ([3]), 
/ij(x, v) and /ij(x, f ) have the same sign (here s), and fi{x) and fi^x+VjCj) are 
on both sides of Xj+fi/2 since fi{x) and fi^x+VjCj) are. In other words, (j, s, i) 
is an edge of Gf{x, f). So (jlj) is proved and, as an immediate consequence, 

UGp)CT,{Gf) {t = l,...,n). (5) 

Then, for all (x, v) G X', we have the following: 

Vertex 1 belongs to none positive circuit of Gp{x,v). (6) 

Indeed, suppose, by contradiction, that vertex 1 belongs to a positive circuit 
of Gp{x, v). Let j be the predecessor of 1 in this circuit, and let t = xi + f i/2. 
By definition, t G Ti{Gp) and from ([5]) it comes that t G Ti{Gf)- We then 
deduce, from ([1]) and the fact that the images of /i are in Y, that fi{x) and 
fi{x + VjCj) are not on both sides of t. In other words, there is no edge from 
j to 1 in Gp{x,v), a contradiction. 



Let A be the set of attractors of T{F) and let 

/ = /\{l}. (7) 

Let (x, v) be any element of X'. Since / is a positive feedback vertex set of 
Gf{x,v) and since Gp{x,v) is a subgraph of Gf{x,v), I is also a positive 
feedback vertex set of Gp{x,v). We then deduce from ([6]) that J is a posi- 
tive feedback vertex set of Gp{x,v). Since this holds for all {x,v) G X', by 
induction hypothesis, 

iAi<n 17^.(^^)1+1, 

i£l 



and from (151) we obtain: 



<n m(G^)i + L 



Now, let A be the set of attractors of T{F), and let Ay be the set of A G A 
containing a point x such that Xi G Y. We claim that: 

VA G Ay, there exists A G A such that A C A. (9) 

So let A G Ay, and consider the set A of x G A such that Xi E Y. We 
prove that A is a trap domain of r(F). Suppose {x, y) to be an edge of r(F) 
such that X G A. By definition, there exists index i such that fi{x) ^ Xi and 
y = X + sign(/i(x) — Xi)ei. We consider two cases: 

(1) Case i > 1. Then, yi = xi & Y. In addition, fi{x) = fi{x) so {x,y) is an 
edge of r(F). Hence y & A (since x & A) and we deduce that y E A. 

(2) Case i = 1. Suppose that Xi < /i(x) (the proof is similar if xi > fi{x)). 
Then, Xi < yi < fi{x) and since Xi and fi{x) are in Y we have t/i G F. 
In addition, min(F) < Xi < /i(x) so Xi < /i(x) < /i(x). Thus (x, y) is 
an egde of T{F). Hence y & A (since x G A) and we deduce that y & A. 

Since y & A in both cases, A is trap domain of T{F). Thus there exists at 
least one attractor A E K such that A C A, and ([9]) holds since A (^ A. 

Following ([9]), there exists a map iJ : Ay — > A such that H{A) C A for all 
A G Ay. Since the attractors of r(F) are mutually disjointed, the elements 
of Ay are mutually disjointed, and we deduces that the images of H are also 
mutually disjointed. Consequently, H is an injection. So |Ay| < |A| and we 
deduce from ([8]) that 

|Ay|<n |T.(Gf)| + L 

is/ 



Since this inequality holds for all y G P, and since A = UygpAy, we have: 

iAi < E i^^i < E [n m{GF)\ + 1] = |P| n m{GF)\ + 1. 

Using ([2]) and ([7]) we conclude: 



|A| < [|ri(G^)| + ij n \T^iGF)\ + 1 = n \T^iGF)\ + 1. 

iei i^i 



u 



Corollary 1 If I is a positive feedback vertex set of G{F), then the number 
of attractors in r(F) and, in particular, the number of fixed points for F are 
less than Yiiei \-^i\- 

Proof — It is sufficient to apply Theorem [5] by noting that: (1) each Gf{x,v) 
is a subgraph of G{F); (2) |Tj(G'f)| + 1 < |Xj|; (3) The number of fixed points 
for F is less than the number of attractors in T{F). D 

Remark 4 The bound on the number of fixed points for F given CorollaryU\ 
has been proved by Aracena et al [IfS] in the Boolean case and under the 
strong hypothesis that G{F) does not contain both a positive and a negative 
edge from one vertex to another (that is, the entries of the Jacobian matrix 
of F are everywhere > or everywhere < 0); see the Theorem\^ stated in the 
introduction. 

Remark 5 Theorems [I] and stated in the introduction are obtained from 
Theorems^ and\^ by noting that Gf{x,v) is a subgraph of Gf{x,v) and by 
using the points (2) and (3) in the proof of Corollary [I]. 



5 Comments 



5. 1 Influence of connections between positive circuits 



Corollary [T] is sufficient to highlight the fact that: "A high level of connection 
between positive circuits leads to a small number of fixed points" . Suppose, for 
sake of simplicity, that all the intervals Xj are of cardinality g, and let r be 
the smallest number of vertices that a positive feedback vertex set of G{F) 
can contain. Then, the smallest upper bound for the number of fixed points 
for F given by Corollary [1] is g'', and the more the positive circuits of G{F) 
are connected, the more r is small. Indeed, let us say that a vertex represents 
a circuit when it belongs to this circuit. Then, r corresponds to the smallest 
number of vertices allowing the representation of each positive circuit. So, 



10 



the more the positive circuit are connected, the more it is possible to choose 
vertices representing a number of positive circuits, and the more r is small. For 
instance, r is always < to the number p of positive circuits that G{F) contains, 
but r < p whenever G{F) has connected positive circuits, and in the extremal 
case where all the positive circuits of G{F) share a same vertex, r = 1. 



5.2 Thomas' logical method 



In practice, the dynamics of a gene network is often modeled from its interac- 
tion graph Q, typically by using the well known Thomas ' logical method [7f8|9] . 
In few words, Thomas associates to ^ a finite state space X and describes the 
behavior of the interactions of Q by logical param,eters . Then, he deduces from 
the value of these parameters a map F from X to itself whose asynchronous 
state transition graph describes a possible dynamics for the network; see [12] 
for a formal presentation. 

This modeling method is coherent with our notion of interaction graph in the 
sense that, for all parameters values, the resulting map F has the property to 
be such that G{F) is a subgraph of Q [13j. So, thanks to Corollary [1], one can 
say, in the total absence of information on the value of the parameters, that 
following Thomas' logical method, the number of attractors in the dynamics 
of the network is less than 

^^{Q,X)= min Hl^^l' 

where I{Q) is the set of smallest positive feedback vertex sets of Q (with 
respect to the inclusion relation). This result is of practical interest since the 
value of the parameters is most often unknown and difficult to estimate, and 
since the number of attractors is an important feature of the dynamics of 
the network. For instance, if the network is known to control a differentiation 
process into k cell types, one often considers that the dynamics of the network 
has to contain at least k attractors. The bound fi{Q,X) can then been used 
in order to check if the data of Q and X is consistent with the presence of k 
attractors (there is inconsistence whenever fi{Q,X) < k). 

5. 3 Feedback circuit functionality 



Finally, Theorem O is related to one of the main concept raised by the Thomas' 
logical method: the concept oi feedback circuit functionality pT|8|9llH] . Roughly 
speaking, it has been observed that some inequality constraints on the logical 
parameters describing the behavior of the interactions of a positive (resp. neg- 
ative) circuit of Q often lead to a dynamics containing several attractors (resp. 



11 



describing oscillations). For that reason, when these constraints are satisfied, 
the corresponding circuit is said functional. Even if this notion is not well un- 
derstand and often informally stated, it is often used in practice to establish 
the value of the logical parameters, see |15fl6fl7lll8fl9f2Uj for instance. 

A natural formalization of the notion of functional circuit, also proposed in 
p!3]l2T] . is the following: given a map F from X to itself whose interaction 
graph G{F) is a subgraph of Q, a circuit C of ^ is functional at {x,v) G X' 
if C is a circuit of Gf{x,v). It is then easy to see that the upper bound for 
the number of attr actors given by Theorem [5] only depends on the localization 
(inside X') and on the connections of the functional positive circuits of the 
system. In our knowledge, this is one of the first mathematical result relating 
the functional circuits of the system to its global dynamical properties (for 
relations between functional circuits and local dynamical properties, see the 
recent parer |21]). 
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